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GEOMETRIC TOMOGRAPHY 
WITH TOPOLOGICAL GUARANTEES 

Omid Amini* Jean-Daniel Boissonnat* Pooran Memari* 



f^ Abstract. — We consider the problem of reconstructing a compact 3-manifold (with boundary) 
CnJ embedded in R'^ from its cross-sections 5 with a given set of cutting planes V having arbitrary 
,_^ orientations. Using the obvious fact that a point x £ "P belongs to the original object if and only 
^ if it belongs to S, we follow a very natural reconstruction strategy: we say that a point a; £ R'^ 
' belongs to the reconstructed object if (at least one of) its n earest point (s) in V belongs to 5. 
^^ This coincides with the algorithm presented by Liu et al. in |LBD^08| . In the present paper, 
^~H we prove that under appropriate sampling conditions, the output of this algorithm preserves the 
homotopy type of the original object. Using the homotopy equivalence, we also show that the 

O reconstructed object is homeomorphic (and isotopic) to the original object. This is the first time 
that 3-dimensional shape reconstruction from cross-sections comes with theoretical guarantees. 

o 



1. Introduction 

Overview. — This paper deals witli tlie reconstruction of 3-dimensional geometric shapes 

^SJ from unorganized planar cross-sections. The need for such reconstructions is a result of the 

*/^ advances in medical imaging technology, especially in ultrasound tomography. In this context, 

the purpose is to construct a 3D model of an organ from a collection of ultrasonic images. 

When the images are provided by free-hand ultrasound devices, the cross-sections of the 

organ belong to planes that are not necessarily parallel. However, it is only very recently 

^3 that reconstruction from unorganized cross-sections has been considered: A very first work 

by Payne and Toga |PT94| was restricted to easy cases of reconstruction that do not require 

branching between sections. In [BG93], Boissonnat and Geiger proposed a Delaunay-based 

algorithm for the case of serial planes, that has been generalized to arbitrarily oriented planes 

in |DP97j and [BM07J. Some more recent work |lJ WC+05|, lLBD+08] can handle the case 

Cd of multilabel sections (multiple materials). Barequet and Vaxman's work |BV09| extends the 

work of fL BD^OS] and can handle the case where the sections are only partially known. 

Most of previous work has been restricted to the case of parallel cross-sections and is mostly 

based on the simple idea of connecting two sections if their orthogonal projections overlap. 
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This paper analyzes a natural generalization of this idea for the case of non parallel sec- 
tions, that has been proposed by Liu et al. in |LBD"'"08j . We prove that under appropriate 
sampling conditions, the connection between the sections provided by this generalization is 
coherent with the connectivity structure of the object and the proposed reconstructed object 
is homeomorphic to the object. To the best of our knowledge, this work is the first to provide 
such a topological study in shape modeling from planar cross-sectional data. The only existing 
results studying the topology of the reconstructed object are restricted to the 2D variant of 



the problem ( |SBG06] and |MB08| ). We study the 2D problem in Section 2.4 of this paper. 
Reconstruction Problem. — Let O C M'^ be a compact 3-manifold with boundary (de- 



noted by dO) of class C"^ ^^' The manifold O is cut by a set V of so-called cutting planes 



that are supposed to be in general position in the sense that none of these cutting planes are 
tangent to dO. For any cutting plane P £ V, we are given the set of intersection O P. 
There is no assumption on the geometry or the topology of these intersections. The goal is to 
reconstruct O from the given set of intersections denoted by S. 

Methodology. — We know that a point x € V belongs to the original object if and only if 
it belongs to S. The goal is now to determine whether a point x G M^ belongs to O or not. 
We follow a very natural reconstruction strategy: 

We say that a point x G M^ belongs to the reconstructed object if (at least one of) 
its nearest point (s) in V belongs to S. 

Different distance function (from V) may be used in order to satisfy properties of interest for 
different applications (for example, to promote the connection between sections in the case of 
sparse data, or to impose a favorite direction to connect the sections, etc). A natural idea is to 
use the Euclidean distance. In this case, the reconstructed object coincides with the method 
introduced by Liu et al. in JLBD"'"08] . We will analyze this method and present appropriate 
sampling conditions providing topological guarantees for the resulting reconstructed object. 

Organization of the paper. — After this brief introduction, in Section [2] we provide 
a detailed description of the reconstructed object TZ. The rest of the paper will be then 
devoted to prove that in the general case, under two appropriate sampling conditions, TZ and 
O are homotopy equivalent, and are more strongly homeomorphic. Indeed, the first sampling 
condition, called the Separation Condition, discussed in Section [2. 1[ ensures good connectivity 
between the sections, but does not necessarily imply the homotopy equivalence. 

As we will see, in order to ensure the homotopy equivalence between TZ and O, a second 



so-called Intersection Condition is required, c.f.. Section 3.4 To make the connection between 



the upcoming sections more clear, in Section 3.1 we shortly outline the general strategy em- 



ployed in proving the homotopy equivalence between TZ and O. In Section [3. 6[ we provide a 
set of properties on the sampling of cutting planes to ensure the Separation and the Inter- 



section Conditions. Finally, in Section 3.8 we show that the two shapes O and TZ are indeed 
homeomorphic (and even isotopic). Some preliminary notions of homotopy theory we use here 
are recalled in Appendix [A) 



'^^This double-smoothness assumption is to simplify a technical part of the homotopy equivalence proof in this 
paper. Indeed, the same topological guarantees hold for a more general case where dO is of class C^'^, i.e., 
when it is continuously differentiable and the normals satisfy a Lipschitz condition. 
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2. Characterization of the Reconstructed Object and First Condition 

As said in the methodology section, we define the reconstructed object 7^ as the set of points 
a; G M'^ so that (at least one of) the nearest point (s) in 7^ to x belongs to S. This is based on 
the distance from the set of cutting planes V. Since we consider the Euclidean distance, this 
involves the arrangement of cutting planes, i.e., the subdivision of M^ into convex polyhedral 
cells induced by the cutting planes. If a point x G M^ belongs to a cell C of this arrangement, 
then its nearest point (s) in V belong(s) to the boundary of C. Thus, we can decompose the 
problem into several subproblems as follows. 

We can restrict our attention to a cell of the arrangement and reduce the reconstruction of 
C to the reconstruction of Oq := O DC for all cells C of the arrangement. 

Reconstructed Object in a Cell of the Arrangement. — We now focus on a cell C of the ar- 
rangement and describe the reconstructed object TZq in C. On each face / of C, the intersection 
of the object O with / is given and consists of a set of connected regions called sections. By 
definition, the sections of a face of C are disjoint. However, two sections (on two neighbor 
faces of C) may intersect along the intersection between their two corresponding faces. The 
boundary of a section A is denoted by dA and is a set of closed curves, called section- contours, 
that may be nested. 

Let us write dC for the boundary of C, and J-^ for the set of faces of C. In the sequel, Sq 
denotes the union of sections of all the faces of C and a point of Sc is called a section-point. 
By the definition of TZ, we have: 

A point x G C is in the reconstructed object TZc if one of its nearest points in dC 
is in Sc. 

The definition of the reconstructed object in a cell C of the arrangement of the cutting 
planes naturally involves the Voronoi diagram of the faces of C defined as follows. 

Voronoi Diagram of a Cell. — For a face / of C, the Voronoi cell of /, denoted by Vorc(/), 
is defined as the set of all points in C that have / as the nearest face of C, i.e., 

Vorc(/) = {xeC\ d{x, f) < d{x, /'), V face /' of C }. 

Here d{., .) stands for the Euclidean distance. The collection of all Vorc(/) when / runs over 
all the faces of C forms a tiling of C, called the Voronoi diagram of C. 





Figure 1. A 2D illustration of the partition of a cell C by the Voronoi Skeleton 
VorSkel(C). Left) The original shape Oc- Right) The reconstructed object TZc- 
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We write dYorcif) for the boundary of Vorc(/). The union of 9Vorc(/) for all the faces / of 
C is called the Voronoi Skeleton of C, and is denoted by VorSkel(C). VorSkel(C) is also called 
the medial axis of the cell, and is the locus of all the points in C that are at the same distance 
from at least two faces of C. To simplify notation, when the cell C is understood from the 
context, we simply remove the index C and write Vor(/), 9Vor(/), etc. 

Definition 1 (Nearest Point Characterization). — For any point x in C, the nearest 
point in dC to x is the orthogonal projection of x onto the nearest face / of C. This projection 
is denoted by npAx). The set of all nearest points to x in dC is denoted by Npq[x). Note that 
for any x ^ VorSkel(C), Npq{x) is reduced to a single point. Based on this, and to simplify 
the presentation, sometimes we drop the index /, and by np{x) we denote a point of Npq{x). 

Definition 2 (Lift Function). — Let a; G C be a point in the Voronoi cell of a face / of C. 
The lift of X in C, denoted by liftc(x) (or simply lift(x) if C is trivially implied), is defined 
to be the unique point of d\orc{f) such that the line defined by the segment [x,lift(x)] is 
orthogonal to /. In other words, lift(a;) is the unique point in 9Vorc(/) whose orthogonal 
projection onto / is np{x). 

— The lift of a set of points X QC, denoted by lift(X), is the set of all the points lift(x) for 
xGX, i.e., lift(X) := { lift(x) \xeX}. 

— The function C : C ^ VorSkel(C) that maps each point a; G C to its lift in VorSkel(C) will 
be called the lift function in the sequel. For any Y C VorSkel(C), C^^iY) denotes the set of 
points X G C such that lift(x) = y for some y £Y. 

We now present a geometric characterization of the reconstructed object using the described 
lifting procedure. 

Characterization of the Reconstructed Object IZq- — If 5c = 0) then, as we said before, for 
any point x G C, np{x) ^ Sq and so IZq is empty. Otherwise, let A G ^c be a section lying on 
a face of C. For each point a ^ A, the locus of all the points x ^ C that have a as their nearest 
point in dC is the line segment [a,lift(a)] joining a to its lift. Therefore, the reconstructed 
object IZc is the union of all the line-segments [a, lift(a)] for a point a in a section A G Sq, 
i.e., 

AeSc aeA 

Note that according to this characterization, if the lifts of two sections intersect in VorSkel(C), 
then these two sections are connected in IZq- This generalizes the classical overlapping criterion 
for the case of parallel cutting planes. The union of all the pieces TZq over all cells C will be 
the overall reconstructed object IZ. This definition is what has been proposed by Liu et al. in 
|LBD+08| . 

The rest of the paper is devoted to prove that under two appropriate sampling conditions, 
TZ and O are homotopy equivalent, and are indeed homeomorphic (and isotopic). 

Let us first infer the following simple observation from the described characterization of 

7^c. 

Proposition 1. — The lift function C : TZq — )• lift(5c) is a homotopy equivalence. 

This is inferred trivially from the fact that the lift function retracts each segment [a, lift(a)] 
onto lift (a) continuously. See Figure p] for an example. 
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Figure 2. A 3D reconstruction example from a pair of parallel sections (in blue). 
The lift function retracts TZc (in green) onto the lift of the sections. 



2.1. First Sampling Condition: Separation Condition. — In this section, we pro- 
vide the first samphng condition, under which the connection between the sections in the 
reconstructed object 7^ will be the same as in the original object O. Our discussion will be 
essentially based on the study of the medial axis, that we define now. 

Definition 3 (Medial Axis of dO) . — Consider dO as a 2-manifold without boundary 
embedded in M^. The medial axis of dO denoted by MA[dO) contains two different parts: 
the so-called internal part denoted by MAi{dO), which lies in O, and the so-called external 
part denoted by MAe{dO), which lies in M^ \ O. 

The internal retract rrii : dO — >• MAj(9C') is defined as follows: for a point x G dO, mi{x) 
is the center of the maximum ball entirely included in O which passes through x. For any 
X G do, mi{x) is unique. Symmetrically, we define the external retract me '■ dO — )• MAe(50): 
for a point x G dO, me{x) is the center of the maximum ball entirely included in M^ \ O which 
passes through x. For any x G dO, me{x) is unique but may be at infinity. In the sequel, we 
may write m{a) for a point in {?7Zj(a),me(o)}. 

The interesting point is that as discussed below if the sample of cutting planes is sufficiently 
dense, then the internal part of MA{dO) lies inside the defined reconstructed object and the 
external part of this medial axis lies outside the reconstructed object. 



Definition 4 (Separation Condition). 

the Separation Condition if 



We say that the set of cutting planes verifies 



MA^{dO)cn and MAeidO) CR^ \n. 

In other words, dTZ separates the internal and the external parts of the medial axis of dO. 
(That is where the name comes from.) 

We will show that in each cell C, the Separation Condition implies that dlZc separates the 
internal and the external parts of the medial axis of dOc- In order to study the Separation 
Condition in a cell C, we will need to consider the medial axis of Oq denoted by MAj(9C'c), 
defined as follows. 
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Definition 5 (Medial axes in a cell C of the arrangement) 

The medial axis of Oq, MAi{dOc), is defined as the set of all points in Oc with at least two 
closest points in dOc, see Figure [3] Symmetrically, by MAf.{dOc) we denote the medial axis 
of the closure of C \ Oq ■ 

Note that the two sets MAi{dOc) and MAi{dO) n C may be different. 



Figure 3. 2D example of medial axes in a cell C of the arrangement. 

We also consider the internal retract rrii^c ■ dOc -^ MAj(OOc) defined as follows. For a 
point X G dOc, f^ifiix) is the center of the maximum ball entirely included in Oc which 
passes through x. Symmetrically, we can define the external retract rriefi : dOc — )• MAe((9C'c): 
for a point x G BOq^ mefi{x) is the center of the maximum ball entirely included in C \ Oq 
which passes through x. It is easy to see that for any x € dO nC, the segments [x,m-j^c(x)] 
and [x^mefi{x)] are subsegments of [x,mi(x)] and [x,me{x)] respectively, and lie on the line 
defined by the normal to dO at x. 

Lemma 1 (Separation Condition Restricted to C). — // the Separation Condition is 
verified, then MAi{dOc) C TZc and MAe{dOc) ^C\nc. 

Proof. — We prove the first part, i.e., MAi{dOc) C TZq- A similar proof gives the second 
part. Let m be a point in MAi{dOc)- Let B{m) be the open ball centered at m which passes 
through the closest points to m in dOc- Two cases can happen: 

• Either, the closest points to m in dOc are in 50, in which case m is a point in MAj((90). 
The Separation Condition states that MAj(c?0) C 7^, and so tti S IZc = TZDC. 

• Otherwise, one of the closest points to m in dOc is a point a in some section A £ Sc- If 
a is on the boundary of A, then since along the section-contours dOc is non-smooth, a lies in 
MAi(dOc) and coincides with m, and m = a is trivially in TZc- Hence, we may assume that a 
lies in the interior of A. Therefore, the ball B{m) is tangent to A at a, and the line segment 
[a, m\ is orthogonal to A. Since B{m) CidC = 9, m and a are in the same Voronoi cell of the 
Voronoi diagram of C. Thus, a G ^c is the nearest point in dC to m. By the definition of TZc, 
we deduce that m G TZc- D 

Assume that the Separation Condition is verified. The first idea which comes to mind is 
to retract points of dO to dTZ by following the normal-directions. A point x £ dO which 
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lies outside 7^ can move towards mi[x) € TZ and stop when dlZ is reached. A point x € dO 
which hes inside TZ, can move toward me{x) and stop when 57?. is reached. According to a 
theorem by Wolter |Wol92j. since dO is assumed to be of class C^, this deformation will 
be a continuous retraction if each normal intersects dTZ at a single point. In such a case, 
do can be deformed to dTZ homeomorphically. But a major problem is that TZ may have a 
complex shape (with cavities), so that a normal to dO intersects dTZ in several points. In such 
a case, such a retraction is not continuous and does not provide a deformation retract of O 
onto TZ. However, we will be essentially following this intuitive idea by looking for a similar 
deformation retract of O onto a subshape of 7^ (the so-called medial shape) . 

2.2. Guarantees on the Connections Between the Sections. — We now show that 
if the sample of cutting planes verifies the Separation Condition, then in each cell C of the 
arrangement, the connection between the sections is the same in Oc and TZc- 

Theorem 1. — If the sample of cutting planes verifies the Separation Condition, TZc o-i^d Oc 
induce the same connectivity components on the sections of C. 

Proof. — The proof is given in two parts: 
(I) If two sections are connected in Oc, then they are connected in TZc. Let A and 
A' be two sections in a same connected component K of Oc ■ Due to the non-smoothness 
of dOc at the boundary of the sections, dA and dA' are contained in MAi{dOc). Thus, 
there is a path 7 in MAi(dOc) D K that connects a point a G dA to a point a' G dA' . 
According to Lemma III MAi{dOc) C TZc. Thus, 7 is a path in TZc that connects A to 
A'. 
(II) If two sections are connected in TZc, then they are connected in Oc. Let A and 
A' be two sections connected in TZc. Let 7 be a path in TZc that connects a point a £ A 
to a point a' G A'. For the sake of contradiction, suppose that a and a' are not in the 
same connected component of Oc. In this case, since 7 joins two points in two different 
connected components of Oc, it must intersect MAe{dOc). But this is a contradiction 
with the fact that 7 C TZc, indeed, according to Lemma [I] we have M.Ae{dOc) riTZc = 0. 

D 

We state the following proposition that will be needed later in this section and in Section [3] 

Proposition 2. — Under the Separation Condition, any connected component of dO is cut 
by at least one cutting plane. 

Proof. — Suppose that X is a connected component of dO which is not cut by any cutting 
plane. There exists a cell C of the arrangement of hyperplanes such that one of the following 
two (symmetric) cases can happen: Either, there exists a connected component H oi O which 
lies in the interior of C such that K C dH, Or, there exists a connected component H of the 
closure of M^ \ O which lies in the interior of C such that K C dH. Without loss of generality, 
suppose that dH bounds a connected component H oi O in C, H is entirely contained in the 
interior of C, and K C dH, see Figure H^ (the other case follows similarly). Take a point m 
in the medial axis of H, i.e., m G H MAi{dO). According to the Separation Condition, m 
belongs to 7^. Thus, by the definition of TZ, one of the nearest points of m in dC, say np{m), 
belongs to S. Since H is not cut by any cutting plane, H D dC is empty and np{m) ^ H. 
Therefore, m and np{m) are in two different connected components of O, and the segment 
[m, np{m)] should intersect MAe{dO) at a point x. On the other hand, by the definition of 
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TZ, the segment [m, np{m)] C TZ. This contradicts the assumption of Separation Condition 
that 7^nMAe(aC') = 0. D 




Figure 4. For the proof of Proposition [2] 



2.3. How to Ensure the Separation Condition? — In this section we provide a sufficient 
condition for ensuring the Separation Condition. For this, we first need some definitions. 

Definition 6 (Reach) . — Let O be a connected compact 3-manifold with smooth boundary 
do in R'^. For a G dO, we define reach(a) = min(d(a, 7raj(a)), (i(a, ?n,e(a))). The quantity 
reach(O) is defined as the minimum distance of dO from the medial axis of dO: 



reachfO) := min d(m,dO) 

m&MA{dO) 



min reach(a). 

a&dO 



Note that since O is compact and dO is of class C^, reach(C') is strictly positive (see |Fed59] 
for a proof). 

Definition 7 (Reach restricted to a cell). — Given a cell C of the arrangement, we de- 
fine reachc(C') = min d{a,m{a)), where either a G dO n C or m{a) G MA{dO) CiC. By 
definition, we have reach(C') = mine (reachc(C')). 

Definition 8 (Height of a Cell). — Let C be a cell of the arrangement of the cutting 
planes. The height of C, denoted by /ic, is defined as the maximum distance of a point x G C to 
its nearest point in the boundary of C, see Figurep] In other words, he := mayixeC d{x, np{x)). 

We remark that the height of any cell C is at most half of the diagonal of C. However, as the 
example of Figure p] (right figure) shows, it may be much smaller than half of the diagonal. 
We now show that by bounding from above the height of the cells by a factor related to the 
reach of the object, we can ensure the Separation Condition. 

Lemma 2 (Sufficient Condition). — // the sample of cutting planes is sufficiently dense 
so that he < reachciO) for any cell C of the arrangement, then the Separation Condition is 
verified. 
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VorSkel(C) 



Diag{Uj' 



w 



3 



np(x) 



Figure 5. Left) Definition of the height of a cell C. Right) Height of a cell C is 
bounded from above by half of the diagonal of C, but may be much smaller for many 
configurations of cutting planes. 



Proof. — The proof is straightforward. Let m be a point in MA(dO) in a cell C of the 
arrangement. We have 

d(?7i, np{m)) < he < reachc(C') < d{m, dO). 

Therefore, the Separation Condition is verified. D 

2.4. Guarantees for 2D Shape Reconstruction from Line Cross-Sections. — We 

proved that under the Separation condition, the connectivity between the sections of C induced 
by the reconstructed object TZq is coherent with the original shape Oq- We presented the proofs 
for the 3D case of the problem, but the proofs easily show that result is valid in any dimension. 
We show in this section that this is strong enough to imply the homotopy equivalence between 
TZq and Oq for the 2-dimensional variant of the reconstruction problem. 

Consider the 2-dimensional variant of the reconstruction problem, that consists of recon- 
structing a 2D-shape from its intersections with arbitrarily oriented cutting lines. In this case 
the sections are line-segments. 

We can focus on a cell C of the arrangement of the plane by the cutting lines. Similar definitions 
for the Voronoi diagram and the Voronoi skeleton of C, the lift function and the reconstructed 
object TZc can be considered. Under the Separation Condition, there is a bijection between 
the connected components of IZc and Oq- By the definition of the reconstructed object, it is 
easy to see that any connected component of IZc is a topological disk. On the other hand, 
according to Proposition [2] under the Separation Condition, any connected component of dO 
is cut by at least one cutting line. We easily deduce that any connected component of Oq is 
a topological disk. 

Using the sufficient condition presented in the last section, if for any cell C of the arrange- 
ment, he < reachc(C'), then the Separation Condition is ensured. We deduce the following 
theorem. 

Theorem 2 (Provably Good 2D Reconstruction) . — // for any cell C of the arrange- 
ment of the cutting lines, he < reachc{0), then IZ is homeomorphic to O. 

Proof. — According to Theorem[T] under the separation condition, there is a bijection between 
the connected components of IZq and Oq. On the other hand, according to the discussion 
preceding the theorem, all the connected components of Oq or IZq are 2-dimensional disks. 
Therefore, there is a homotopy equivalence between each pair of corresponding connected 
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components of Oq or 7^^. This provides a homotopy equivalence between TZc and Oc- As we 



will explain in detail in Section 3.5, the homotopy equivalences in the different cells of the 
arrangement can be extended to a homotopy equivalence between TZ and O. Finally, since TZ 
and O are two homotopy equivalent 2-dimensional submanifolds of E?, we deduce that there 
is a homeomorphism between 7^ and O. D 



3. Topological Guarantees 

To clarify the connection between the upcoming sections, let us shortly outline the general 
strategy employed in proving the homotopy equivalence between 7^ and O. 



3.1. Proof Outline of the Homotopy Equivalence Between TZ and O. — We will 
provide a homotopy equivalence between TZc ^i^d Oc in each cell of the arrangement. (And then 
glue these homotopy equivalences together to form a global homotopy equivalence between 7^ 



and O.) In Section 2.2 we showed that under the first sampling condition called Separation 
Condition the connection between the sections in the reconstructed object TZc is the same as 
in Oc, in the sense that there is a bijection between the connected components of TZc and 
the connected components of Oc- This implies that for proving the homotopy equivalence 
between TZc and Cci it will be enough to show that the corresponding connected components 
have the same homotopy type. In order to extend these homotopy equivalences to a homotopy 
equivalence between TZ and O, we will have to glue together the homotopy equivalences we 
obtain in the cells of the arrangement. This needs some care since the restriction to a section 
S of the two homotopy equivalences defined in the two adjacent cells of S may be different. 
To overcome this problem, we need to define an intermediate shape A4c in each cell C, called 
the medial shape. The medial shape has the following three properties: 

(i) The medial shape contains the sections of C, i.e., Sc C A4c- 

(ii) There is a (strong) deformation retract r from Oc to Aic- In particular, this map is 
a homotopy equivalence between Oc and A4c- And its restriction to Sc is the identity 
map. 
{iii) Under the first sampling condition (Separation Condition), Aic ^ ^C- 

The first two properties will be crucial to guarantee that the homotopy equivalences conform 
on each section under the Separation Condition. Indeed, the map Oc — ?■ A^c ^^ TZc, obtained 
by composing the deformation retract and the inclusion, restricts to the identity map on each 
section of Sc- Thus, we can glue all these maps to obtain a global map from O to TZ. 



Using a generalized version of the nerve theorem (see Section 3.5) and property {ii) above, 
we can then reduce the problem to proving that the inclusion i : Mc ^^ T^C forms a homotopy 
equivalence in each cell. Using Whitehead's theorem, it will be enough to show that the 
inclusion i induces isomorphisms between the corresponding homotopy groups. Under the 
Separation Condition, we prove that i induces an injective map on the first homotopy groups, 
and that all higher dimensional homotopy groups of Mc and TZc are trivial. Unfortunately, 
the Separation Condition does not ensure in general the surjectivity of i on the first homotopy 
groups. To overcome this problem, we need to impose a second condition called Intersection 
Condition. Under the Intersection Condition, the map i will be surjective on the first homotopy 
groups, leading to a homotopy equivalence between O and 7^. 
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According to the guarantees on the connectivity between the sections (Theoremnj), to prove the 
homotopy equivalence between TZq and Oq under the Separation Condition, we may restrict 
our attention to each of the corresponding connected components. 

In the sequel, to simplify the notations and the presentation, we suppose that Oq 
and thus TZc are connected, and we show that Oq and TZc have the same homotopy 
type. It is clear that the same proofs can be applied to each corresponding connected 
components of Oq and IZq to imply the homotopy equivalence in the general case of 
multiple connected components. 

3.2. Medial Shape. — We now define an intermediate shape in each cell C of the ar- 
rangement called the medial shape. The medial shape enjoys a certain number of important 
properties, discussed in this section, which makes it playing an important role in obtaining 
the homotopy equivalence between Oc and TZq- 

Definition 9 (Medial Shape M.c)- — Let a; be a point in Sc C dOc- Let w{x) = 
[x,mifi{x)\ be the segment in the direction of the normal to dOc at x which connects x to 
the point mifi{x) G MAj (c^Oc)- We add to MAj(9C'c) all the segments w{x) for all the points 
x ^ Sc- We call the resulting shape A^Ci see Figure p] (Left) for an example. More precisely, 

Mc:=MAi{dOc)U{\J w{x)). 
xeSc 




Figure 6. A 2D illustration of the medial shape Aic in purple (left) and Aic in red (right). 

Proposition 3. — The medial shape verifies the following set of properties: 

(i) The medial shape contains the sections of C, i.e., Sq ^ A^c- 

{a) There is a (strong) deformation retract r from Oc to Mq- In particular, this map is a 
homotopy equivalence between Oc and Mq- ^nd its restriction to Sq is the identity map. 
{Hi) Under the Separation Condition, Mq ^ ^C- 

Proof. — (i) This property is true by the definition of the medial shape. 

(a) This is obtained by deforming Oq to Aic in the direction of the normals to the boundary 

dOc. Note that the boundary dOc is smooth except on the boundaries of sections in 

Sc, and the boundaries of the sections in Sc are already in Aic, thus, the deformation 

retract is well defined. Moreover, since dO (and so dOc) is supposed to be of class C^, 
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according to a theorem by Wolter |Wol92| . this defo rmation is continuous and is hence 
a continuous deformation retract from Oc to A^c 



(2) 



{iii) Since Mc = MAi{dOc) U (U^g-Sc '^(^)) ^"^ ^" addition MAi{dOc) C 7^c, it will be 
sufficient to show that for any x in a section A £ Sq, w{x) C TZq- (Recall that w{x) 
is the orthogonal segment to dOc at x that joins x to the corresponding medial point 
mifi{x) in MAi{dOc)-) We will show that w{x) is contained in the segment [x,lift(a;)]. 
The point x is the closest point in dOc to mifi{x). Thus, the ball centered at mifi{a) 
and passing through x is entirely contained in O and its interior is empty of points 
of dC. Thus, in the Voronoi diagram of C, mic{a) is in the same Voronoi cell as x. 
On the other hand, x is the closest point in Sq C dOc to lift(x). It easily follows 
that d{x,\\it{x)) > d{x,mifi{x)). It follows that the segment [x,mifi{x)] = w{x) is a 
subsegment of [a;, lift(x)]. Therefore, by the definition oiTZc, w{x) C TZc- 

D 

We end this section with the following important remark and proposition which will be 
used in the next section. By replacing the shape Oc with its complementary set we may 
define an exterior medial shape Mq- This is more precisely defined as follows. Let O be the 
closure of the complementary set of O in M^. And let Oc be the intersection of O with the 
cell C. The medial shape of Oq-, denoted by Mqi is the union of the medial shapes of the 
connected components of Oc, see Figurep (Right). Similarly, under the Separation Condition, 
the following proposition holds. 

Proposition 4- — Let Oq be the closure of the complementary set of Oq in C and Mq be 
the medial shape of Oc- Under the Separation Condition: (i) There is a strong deformation 
retract from C \ Aic ^o Oq, and [ii) We have IZq C C \ Mc- 

Proof. — The proof of Property (i) is similar to the proof of Proposition ^ by deforming 
along the normal vectors to the boundary of Oq- The second property {ii) is equivalent to 
A^c C C \ 7^c. D 

3.3. Topological Guarantees Implied by the Separation Condition. — Throughout 
this section, we suppose that the Separation Condition holds. By the discussion at the end 



of Section 3.1 and without loss of generality, we may assume that Oc and hence TZc are 
connected. Thus, Oc and IZc are connected compact topological 3-manifolds embedded in M^. 

In Section [3^ we defined the medial shape Mc and showed that Mc is homotopy equivalent 
to Oc by giving a (strong) deformation retract from Oc to Mc ■ We have also shown that under 
the Separation Condition, Mc C TZc- Using these properties, the goal will be to prove that the 
inclusion i : Mc ^^ TZc is a homotopy equivalence. As the objects we are manipulating are all 
CW-complexes, homotopy equivalence is equivalent to weak homotopy equivalence according 
to Whitehead's theorem. 

Hence, it will be enough to show that i : Mc ^^ TZc induces isomorphism between the 
corresponding homotopy groups. 



"^Indeed, it is possible to show that the same result holds for a more general case where dO is of class C^'^ 
(i.e., it is continuously differentiable and its normal satisfies a Lipschitz condition). In that case, Welter's 
theorem on the continuity is not valid anymore (i.e., there are counter examples). However, one can prove the 
existence of a continuous deformation in our spacial case by using Lieutier's deformation based on generalized 
distance function. We refer the interested reader to Lieutier's paper |Lie04| for more details. 
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Injectivity on the Level of Homotopy Groups. — We will make use of the lift function in C, 
c.f. Section [2] We consider the restriction of the lift function to the reconstructed object TZq- 
According to Proposition 111 the lift function C : TZc -^ lift (5c) is a homotopy equiva- 
lence. We claim that it will be sufficient to show that the restriction of the lift function to 
Aic is a weak homotopy equivalence. Because in this case, using the following diagram and 
Whitehead's theorem, we can infer that i : Aic "^ T^c is a homotopy equivalence. 



> TZc 




lift (5c) 

More precisely, if C : Mq ~^ lift (5c) is a weak homotopy equivalence (what we will prove 
below), since C : TZc ~^ lift(5c) is also a homotopy equivalence and because of the commuta- 
tivity of the above diagram, the inclusion i : Mc ^^ TZc induces isomorphisms between the 
homotopy groups of Mc and TZc- Thus, Whitehead's theorem implies that i is a homotopy 
equivalence. We first show that under the Separation Condition, the restricted lift function 
C\^^ induces injections on the level of homotopy groups. 

Theorem 3 (Injectivity) . — Under the Separation Condition, the homomorphisms between 
the homotopy groups of Aic f^-i^d lift(5c), induced by the lift function C, are injective. 

Proof. — Under the Separation Condition, we have A4c C TZc. Let Mc be the medial shape 
of the closure of the complementary set of Oc in C. We refer to the discussion at the end of 
the previous section for more details. Recall that by Proposition kl we have TZc C C \ Adci 
and there exists a deformation retract from C \ Mc to Oc (in particular Oc and C \ Aic 
are homotopy equivalent). We have now the following commutative diagram in which every 
map (except the lift function C) is an injection (or an isomorphism) on the level of homotopy 
groups. 

Oc 



Mc 




> TZc ^ > C\Mc 



lift (5c) 



Using this diagram, the injectivity on the level of homotopy groups is clear: For any integer 
J > 1, consider the induced homomorphism C^: : vrj(A^c) ~^ "^j (lift (5c)). Let x S Trj(A4c) 
be so that >C*(a;) is the identity element of 7rj(lift (5c)). It is sufficient to show that x is the 
identity element of TTj^Aic)- Following the maps of the diagram, and using the homotopy 
equivalence between lift (5c) and TZc, we have that z^,(x) is mapped to the identity element 
of -KjiTZc). Then, by the inclusion TZc ^~^ C\ Mc, it goes to the identity element of C \ M.c, 
and by the two retractions, it will be mapped to the identity element of ■7rj{Aic)- As this 
diagram is commutative, we infer that x is the identity element of Trj{Mc)- Thus, £* : 
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TTj{Mc) ~^ '''"j(lift(6'c)) is injective for all j > 1. The injectivity for j = is already proved in 
Theorem [T] D 

We have shown that under the Separation Condition, the lift function C : Aic ~^ lift (5c) 
induces injective niorphisms between the honiotopy groups of A4c and lift(6'c). If these in- 
duced morphisms were surjective, then C would be a honiotopy equivalence (by Whitehead's 
theorem) . We will show below that the Separation Condition implies the surjectivity for all the 
homotopy groups except for dimension one (fundamental groups). Indeed, we will show that 
under the Separation Condition, all the i-dimensional homotopy groups of Mc and lift (5c) 
for i > 2 are trivial. Once this is proved, it will be sufficient to study the surjectivity of 
£*:^i(A^c)^^i(lift(5c)). 

Remark 1. — Note that the injectivity in the general form above remains valid for the cor- 
responding reconstruction problems in dimensions greater than three. However, the vanishing 
results on higher homotopy groups of Oq and IZq are only valid in dimensions two and three. 

The topological structures of TZq and Oq are determined by their fundamental groups. — In 
this section, we show that if the Separation Condition is verified, then the topological structure 
of the portion of O in a cell C (i.e., Oq) is simple enough, in the sense that for all i > 2, the 
i-dimensional homotopy group of Oc is trivial. We can easily show that TZc has the same 



property ^^' As a consequence, the topological structures of Oq and T^c are determined by 
their fundamental group, 7ri(C'c) and 7ri(7^c)- 

We first state the following general theorem for an arbitrary embedded 3-manifold with 
connected boundary. 

Theorem 4- — Let K be a connected 3-manifold in M^ with a (non-empty) connected bound- 
ary. Then for all i >2, iTi{K) = {0}. 

This theorem can be easily obtained from Sphere Theorem, see e.g.. Corollary 3.9 of |Hat02] 
or for more details |MemlO] . From this theorem, we infer the two following theorems. 

Theorem 5. — Under the Separation Condition, TTi{Oc) = {0}, for all i > 2. 

Proof. — We only make use of the fact that under the Separation Condition, any connected 
component of dO is cut by at least one cutting plane (Proposition p|. In this case, every 
connected component of Oq is a 3-manifold with connected boundary. The theorem follows 
as a corollary of Theorem [4] D 

Theorem 6. — 'Ki{lZc) = {0}, for all i >2. 

Proof. — Using Theorem |4j it will be sufficient to show that the boundary of any connected 
component K of TZc is connected. Let x and y be two points on the boundary of K, and let 
S and S' be two sections so that x £ [a, lift(a)] for some a £ S and y £ [6, lift(6)] for some 
b £ S'. By the definition of TZc, x is connected to S in &TZc, and y is connected to S' in dTZc. 
On the other hand, since S and S' are two sections in the connected component K of TZc, 
they lie on dK and are connected to each other in dK (and so in dTZc). Thus, x is connected 
to y in dTZc. D 



"' Recall that for simplifying the presentation, we assume that Oc and so TZc are connected. The same proof 
shows that in the general case, the same property holds for each connected component of Oc or TZc- 
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3.4. Second Condition: Intersection Condition. — In the previous section, we saw 
that under the Separation Condition, the topological structures of Oq and TZq are determined 
by their fundamental group 7ri(Oc) and Tri{TZc), respectively. The goal of this section is to 
find a way to ensure an isomorphism between the fundamental groups of TZq and Oq- We 
recall that as Oq and Aic a-re homotopy equivalent, Tri{Oc) is isomorphic to iTi{Aic)- On the 
other hand, TZq and lift(5c) are homotopy equivalent, and iriiTZc) is isomorphic to 7ri(lift(5c)) 
(c.f. last diagram). Thus, it will be sufficient to compare TTi{Aic) and 7ri(lift(5c)). 

We consider C^, : vri(7Wc) — )■ 7ri(lift(5c)), the map induced by the lift function from A4c to 
lift (5c) on fundamental groups . We showed that £* is injective. A natural condition to ensure 
that C^, is an isomorphism is to impose that lift (5c) is contractihle (or more generally, each 
connected component of lift(5c) is contractible) . This is very common in practice, where the 
sections are contractible and sufficiently close to each other. In this case, all the homotopy 
groups 7rj(lift(5c)) are trivial and by injectivity of £* proved in the previous section, £* 
becomes an isomorphism in each dimension. Hence, the homotopy equivalence between TZq 
and Oc can be deduced. 




Figure 7. Left) Original shape. Right) Reconstructed shape. 3D example of the 
case where the Hft function from M.c to lift(5c) fails to be surjective: xi and X2 are 
two points with the same lift in lift (5c). The lift of any curve 7 connecting xi and 
X2 in AAc provides a non-zero element of 7ri(lift(5c),a;). The reconstructed shape is 
a torus and is not homotopy equivalent to the original shape which is a twisted 
cylinder. 



However, the map C-t : '/ri(7Wc) ~^ tti (lift (5c)) fails to be surjective in general (where the 
connected components of lift (5c) are not necessarily contractible). Figure^ shows two shapes 
with different topologies, a torus and a {twisted) cylinder, that have the same (inter) sections 
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with a set of (two) cutting planes. Hence, whatever is the reconstructed object from these 
sections, it would not be topologically consistent for at least one of these objects. In particular, 
the proposed reconstructed object (7^) is a torus which is not homotopy equivalent to the 
(twisted) cylinder (O). In addition, we note that the Separation Condition may be verified 
for such a situation. Indeed, such a situation is exactly the case when the injective morphism 
between the fundamental groups of O and TZ is not surjective. This situation can be explained 
as follows: let xi and X2 be two points in the sections Si and 5*2 with the same lift x in 
lift (5c)- The lift of any curve 7 connecting xi and X2 in Mc provides a non-identity element 
of 7ri(lift(5c),x) which is not in the image of C^,. We may avoid this situation with the 
following condition. 

Definition 10 (Intersection Condition). — We say that the set of cutting planes verifies 
the Intersection Condition if for any pair of sections Si and Sj in 5c, and for any connected 
component X of lift(S'j) Pi lift(S'j), the following holds: there is a path 7 C Mq from a point 
a £ 5j to a point b G Sj with lift(a) = lift(6) = x G X so that i2*(7) is the identity element of 
VTi (lift (5c ), x), i.e., is contractible in lift(5c) with a homotopy respecting the base point x. 



In Section [3. 6[ we will show how to verify the Intersection Condition. Let us first prove the 
surjectivity of the map d, which is deduced directly from the Intersection Condition. 

Theorem 7. — Under the Intersection Condition, the induced map £* : iri{A4c) ~^ 
7ri(lift (5c)) is surjective. 

Proof. — Let yo be a fixed point of A4c and xq = C{yo). We show that C^ : vri(7Wc)2/o) ~^ 
TTi (lift(5c), Xq) is surjective. Let a be a closed curve in lift(5c) which represents an element of 
TTi (lift (5c ), Xq). We show the existence of an element /3 G '^i{-^Cjyo) such that >C*(/3) = [a], 
where [a] denotes the homotopy class of a in vri (lift (5c ), xq). See Figure 8[ left. 
We can divide a into subcurves ai, . . . , am such that aj joins two points Xj-i and Xj, and 
is entirely in the lift of one of the sections Sj, for j = 1, . . . , m. We may assume yo G Si = 
Sm- For each j = 1,. . . ,m, let /3j be the curve in Sj joining two points Zj to Wj which 
is mapped to aj under C Note that Wj and Zj^i (possibly) live in two different sections, 
but have the same image (xj) under the lift map C Let Xj be the connected component 
of lift(S'j) n lift(S'j-i-i) which contains Xj, see Figure fright. According to the Intersection 
Condition, there is a path 7^ C Mc connecting a point aj £ Sj to a point bj+i G Sj+i such 
that lift(aj) = lift(6j+i) = x'- G Xj and the image of 7^ under C is the identity element of 
VTi (lift (5c ), x') (i.e., is contractible with a homotopy respecting the base point x'). Since Xj 
is connected, there is a path from Xj to x'- in Xj, so lifting back this path to two paths from 
Wj to aj in Sj and from 6j+i to Zj^i and taking the union of these two paths with 7^, we infer 
the existence of a path 7' C Aic connecting Wj to Zj^i, such that the image of 7' under C is 
contractible in lift (5c) with a homotopy respecting the base point Xj. 

Let /3 be the path from xq to xq obtained by concatenating f3j and 7' alternatively, i.e., 
P = /3i7^/3272 • • • /3m-i7m/^m7m- We claim that >C*([/3]) = [a]. This is now easy to show: 
we have JC^{f3) = ai£*(7^)a2 • • ■ ^*{'y'm)'^m, and all the paths vC*(7j) are contractible to the 
constant path [xj] by a homotopy fixing Xj all the time. We deduce that under a homotopy 
fixing Xq, ai£*(7^) . . . £*(7^)am is homotopic to aiQ2 • ■ • am = a, and this is exactly saying 
that £*([/3]) = [a]. The surjectivity follows. D 

Putting together all the materials we have obtained, we infer the main theorem of this section. 
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Figure 8. For the proof of Theorem l7| 3D illustration of two sections Sj and Sj+i 
(in blue) and their lift (in red). For any closed curve a in lift(iSc), we show that there 
exists a closed curve (3 in A4c so that the lift of /3 is honiotopic to a in lift (5c). Left) 
We show the existence of /3 such that Ci,{ji) differs from a only in some contractible 



Right) Notations used in the proof: We divide a into subcurves {aj}™ such 
; is entirely in the lift of the section Sj . By Xj we denote the point shared by 



paths, 
that a 
aj and ctj+i which is the lift of two points Wj G Sj and Zj+i G S'j+i. 



Theorem 8 (Main Theorem-Part I). — Under the Separation and the Intersection Con- 
ditions, IZc is homotopy equivalent to Oq for any cell C of the arrangement. 

3.5. Generalized Nerve Theorem and Homotopy Equivalence of TZ and O. — In 

this section, we extend the homotopy equivalence between TZc and Oc, in each cell C, to a 
global homotopy equivalence between TZ and O. To this end, we make use of a generalization 
of the nerve theorem. This is a folklore theorem and has been observed and used by different 
authors. For a modern proof of a still more general result, we refer to Segal's paper |Seg68| . 
(See also [May 03] , for a survey of similar results.) 

Theorem 9 (Generalized Nerve Theorem). — Let H : X ^ Y be a continuous map. 
Suppose that Y has an open cover /C with the following two properties: 

— Finite intersections of sets in /C are in fC. 

— For each U £ IC, the restriction H : H~^(U) -^ U is a weak homotopy equivalence. 

Then H is a weak homotopy equivalence. 

Let Hq : Oc -^ TZc be the homotopy equivalence obtained in the previous sections between 
Oc and TZc- (So He is the composition of the retraction Oq — >• Mc and the inclusion Mq ^^ 
TZc-) Let H : O ^ TZ he the map defined by H{x) = Hc{x) if x G Oq for a cell C of the 
arrangement of the cutting planes. Note that H is well-defined since Hc\sc — 'i^Sci foi' ^ C. 
In addition, since for all cell C, Hq is continuous, H is continuous as well. 
We can now apply the generalized nerve theorem by the following simple trick. Let e be an 
infinitesimal positive value. For any cell C of the arrangement of the cutting planes, we define 
Oq = {x G M^, d(x, Oc) < e }■ Let us now consider the open covering /C of O by these open 
sets and all their finite intersections. It is straightforward to check that for e small enough, the 
restriction of H to each element of /C is a weak homotopy equivalence. Therefore, according 
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to the generalized nerve theorem, H is a weak homotopy equivalence between TZ and O. And 
by Whitehead's theorem, H is a homotopy equivalence between TZ and O. Thus, we have 

Theorem 10 (Main Theorem- Part II). — Under the Separation and Intersection Condi- 
tions, the reconstructed object TZ is homotopy equivalent to the unknown original shape O. 



3.6. How to Ensure the Intersection Condition? — In Section 2.2 we showed that the 
Separation Condition can be ensured with a sufficiently dense sample of cutting planes. In 
this section we provide a sufficient condition that implies the Intersection Condition. 

We showed that by bounding from above the height of the cells by the reach of the object, we 
can ensure the Separation Condition. In order to ensure the Intersection Condition, we need 
a stronger condition on the height of the cells. As we will see, this condition is a transversality 
condition on the cutting planes that can be measured by the angle between the cutting planes 
and the normal to dO at contour-points. 

Definition 11 (Angle a^). — Let a be a point on the boundary of a section A ^ Sq on 
the plane Pa- We consider mi{a), that may be outside the cell C. We define Oa as the angle 
between Pa and the normal to dO at a, i.e., q^ := angle{PA, [a, rn-j(a)]), see Figure^ 




Figure 9. Definition of a^ (a 2D illustration). 

Sufficient Conditions. — We now define the sampling conditions on the cutting planes. (See 
Section [2. 3| for the definitions of he and reachc(O).) 

(CI) Density Condition : For any cell C of the arrangement, he < reachc(C'). 
(C2) Trans versahty Condition : For any cell C, 

he < - [l — sin(aa)) reach(a), Va G dSe- 

The Density Condition is based on the density of the sections. The Transversality Condition is 
defined in a way that the transversality of the cutting planes to dO and the distance between 
the sections are controlled simultaneously. (Indeed, the significance of sin(aa) is to control 
the transversality, and bounding from above he allows us to control the distance between the 
sections.) 
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Remark on the Transversality Condition. — The transversality of the cutting planes to 
do seems to be a reasonable condition in practice, specially for applications in 3D ultrasound. 
Indeed, according to [Rou03,, Section 1.2.1], from a technical point of view if a cut is not 
sufficiently transversal to the organ, the quality of the resulting 2D ultrasonic image is not 
acceptable for diagnosis. 

According to Lemma [2] of Section |2.2[ the Density Condition implies the Separation Con- 
dition. We will show that under the Transversality Condition, the Intersection Condition is 
verified. Therefore, by increasing the density of the sections of O, with preferably transversal 
cutting planes, we can ensure the required sampling conditions, and as a consequence, provide 
a topologically consistent reconstruction of O. 

Theorem 11 (Main Theorem- Part III). — // the set of cutting planes verifies the Den- 
sity and the Transversality Conditions, then the Separation and the Intersection Conditions 
are verified. Therefore, the proposed reconstructed object TZ is homotopy equivalent to the 
unknown original shape O. 



3.7. Proof of Theorem \11\ — This section is devoted to the proof of Theorem 11 It 



remains to show that under the Transversality Condition, the Intersection Condition is verified. 
We need the following notations. 

Notation i^Ki{Sc) and Ke{Sc)). — Recall that VorSkel(C) is the locus of the points with 
more than one nearest point in dC. We write Ki{Sc) (resp. Ke{Sc)) for the set of points 
X G VorSkel(C) such that all the nearest points of x in dC lie inside (resp. outside) the 
sections. 

Notation {fni{a) and rne{a)^' — Let a be a point on the boundary of a section A ^ Sq 
on the plane Pa- We write rni[a) (resp. rne{a)) for the orthogonal projection of mi{a) (resp. 
me{a)) onto Pa- See Figure 10 Note that we have d{mi{a),fni(a)) = sin(aa) d{a,mi{a)) and 
d{me{a),rne{a)) = sin(aa) d{a,me{a))- 




FiGURE 10. Definition of mj(a) (a 2D illustration). 



Lemma 3- — If the Transversality Condition is verified, then for any a £ dSc, we have 
lift(mj(a)) G Ki{Sc) and lift (me (a)) G Ke{Sc)- In addition, the two segments that join 
lift(mi(a)) to its nearest points in dC both lie entirely in Mq- 
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Proof. — We first show that hft(m.j(a)) G Ki{Sc)- By symmetry, the analogue property for 
lift (me (a)) can be proved similarly. To simplify the notation in the proof of this lemma, we 
write rrii for mi{a), and mi for rni{a), and simply write B(mi) for the ball centered at rrij 
of radius d{mi,a). We have d{mi,rni) < d{mi,a). Thus, rUi lies in B{mi). As this ball is 
contained in O, mj is in C Consider now lift(mj) on VorSkel(C), and call y the point distinct 



from nii such that lift(mi) = lift(y), see Figure 10 To have lift(mj) € Ki{Sc), we need to 
show that y is in O. We have 

d{mi,y) < d{mi,rni) + d{mi,\ift{mi)) + d{lih{y),y) 
< sm{aa) d{a,'mi) + 2 he < d{a,mi). 

Thus, y belongs to B{mi), and we deduce that y £ O and lift(mi) G Ki(Sc)- In addition, the 
ball B centered at lift(mj) which passes through rfii and y is entirely contained in B{mi) C O. 
Thus, the interior of B is empty of points of dOc and B is & medial ball of Oq, and its center 
lift(mj) belongs to MAi(9C'c)- Since rrli and y are in 5c, according to the definition of A^Ci 
the line-segments [lift(mi), mj] and [lift(mj),y] lie entirely in Mq- D 

Let Sa and Sh be two sections in Sq such that lift(S'a) n lift(S';,) is non-empty. (As it 
will become clear, a and b will be two chosen points in Sa and Sb respectively.) We recall 
that C-t denotes the homomorphism £* : ttj{M.c) ~^ ^j(lift(5c)) induced by the lift function 
C : Mc ~^ lift(5c). We will show that for any two points a E dSa,b G dSb such that 
lift(a) = lift(6), there exists a path 7 C A4c between a and b so that i2*(7) is contractible in 
£*(7ri(A^c))- Let Pa and Pb be the cutting-planes of Sa and 5;, respectively. One of the two 
following cases happen: 

• If the segment [a,mi(a)] (i.e., the projection of the segment [a,mi{a)] onto Pa) is not cut 
by any other cutting plane, then we claim that lift (a) is connected to lift(mj(a)) in Ki{Sc)- 
Consider the two-dimensional ball B in the plane Pa centered at rni{a) and passing through 
a. This ball is contained in the 3D ball centered at mi{a) and passing through a, which lies 
entirely inside O. Considering the intersection of this 3D ball with Pq, we infer that B lies 



in the section Sa- Therefore, lift(B) is entirely contained in lift(S'a). In Figure 11 lift(i?) is 
colored in green. 




Figure 1 1 . For the proof of Theorem [TT] the lift of two sections Sa and Sb on 
VorSkel(C) is illustrated. We prove that the lift of the segment [a.mi{a)] lies in 

lift(5a) n lift(56). 
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Since lift(a) = lift(6), we have lift(a) G lift(S'a) n lift(S';,). On the other hand, according to 
Lemma l3| lift(mj(a)) hes in hft(S'a) PI hft(S';,). For the sake of a contradiction, suppose that 
hft(a) is not connected to lift(mj(a)) in Ki(Sc)- In this case, lih{B) intersects (at least) two 



different connected components of Uft(S'a) n hft(S'f,), see Figure 11 Since Hft(-B) C hft(5a). 



we deduce that lift(i?) intersects Hft(5'ft) in (at least) two different connected components. 
Consider now the maximal open ball contained in lift(i?) which is empty of points of lift(5'ft). 
Such a ball consists of the lift of a 2D ball B' in the plane P^ which is tangent to dS^ at two 
points X and x'. B' is the intersection of the 3D medial ball of the complementary set of O 
passing through x and x' . Thus, rne(x), which is the projection of the center of this 3D ball 
(me(x)) onto Pfe, lies in B' . We have lift(me(x)) G lift(P') C Mt{B) C lift(5a). Thus, one of 
the nearest points of lift (me (x)) in dC lies in Sa C Sc- This contradicts lift(me(a)) G Ke(Sc) 
(Lemma ^. 

Therefore, lift (a) is connected to lift(mj(a)) in Ki{Sc)- Let us call a' and b' the nearest points 
of lift(mj(a)) in Sa and Sb respectively, see Figure 12 left. According to Lemma pj the line- 
segments [o',lift(mj(a))] and [5', lift (mj (a))] lie inside Mq- We now define the path 7, colored 
in red in Figure [l2]-left, as the concatenation of four line-segments: [a, a'] C Sa, [a' , lift(mj(a))], 
[6', lift (mj (a))] and [6', 6] C Sh- We know that [a', lift (m^ (a))] and [6', lift (mj (a))] are mapped 
to a single point lift(mi(a)) by the lift function. Thus, the image of 7 under the lift function 
is the line-segment [lift (a), lift (m^ (a))], which is trivially contractible in £*(7ri(A^c)). 





Figure 12. A 2D illustration of the two cases considered in the proof of Theorem [TTj 
the path between a and h in Mq is given in red. 



• If the segment [a, 7raj(o)] (i.e., the projection of the segment [a, mj(a)] onto Pa) is cut by 
a cutting plane P (Figure 12 (Right)), then we claim that a is connected to b in Sq- Let I 



be the lift of the segment [a,mj(a)]. For any point j; G I, we write B{x) for the open ball 
centered at x and tangent to the planes of sections Sa and S^. Since lift(a) G lift(S'(i)nlift(5;,), 
i?(lift(a)) is empty of points of dC. In particular, it does not intersect P. Hence, there exists 



a point xq G I, such that P is tangent to B{xq). See Figure 12 (Right). We also observe that 



P intersects [a,mi{a)\ and the corresponding segment in Sh which is also mapped to / by the 
lift function (the segment [6, b'] in Figure 



12 



(Right)). In this case, since Sq C Mq, there is a 
path in Sc (and so in M.c) between a and 6, colored in red in Figure [l2] (Right) , such that its 
image by the lift function is the concatenation of [lift (a), xq] with a contractible path from xq 
to Xq. 



This ends the proof of Theorem 1 1 
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3.8. Deforming the Homotopy Equivalence to a Homeomorphism. — Using the 
honiotopy equivalence between TZ and O, one can sliow tlie following theorem. 

Theorem 12 (Main Theorem- Part IV). — Under the Separation and the Intersection 
Conditions, the two topological manifolds IZ and O are homeomorphic (in addition, they are 
isotopic). 

Although this result is stronger than the honiotopy equivalence, the way our proof works 
makes essentially use of the topological study of the previous sections. 

Proof. — Again, we first argue in each cell of the arrangement and show the existence of a 
homeomorphism between Oc and TZc whose restriction to Sc is the identity map. Gluing these 
homeomorphisms together, one obtains a global homeomorphism between TZ and O. Let C be 
a cell of the arrangement of the cutting planes. A similar method used to prove the homotopy 
equivalence between TZc and Oc shows that dTZ H C and dO D C are homotopy equivalent 
2-manifolds and are therefore homeomorphic, and in addition there exists a homeomorphism 
f3c ■ do n C — )• dTZ n C which induces identity on the boundary of sections in Sc- We showed 
that the topology of TZc and Oc is completely determined by their fundamental groups, i.e., 
all the higher homotopy groups of TZc and Oc are trivial. Moreover, there is an isomorphism 
between iti{Oc) and 7ri{TZc), and the induced map (/3c)* : TTi{dO Pi C) — t- 'Ki{dlZ Pi C) on 
first homotopy groups is consistent with this isomorphism (in the sense that there exists a 
commutative diagram of first homotopy groups). This shows that there is no obstruction 
in extending (3c to a map ac '■ Oc — ^ TZc inducing the corresponding isomorphism between 
iTi{Oc) and 7ri(7^c)) and such that the restriction of ac to Sc remains identity. Since all 
the higher homotopy groups of Oc and TZc are trivial, it follows that ac is a homotopy 
equivalence. We can now apply the following theorem due to Waldhausen which shows that 
a can be deformed to homeomorphism between Oc and TZc by a deformation which does not 
change the homeomorphism ac between the boundaries. A compact 3-manifold M is called 
irreducible if 7r2(M) is trivial. We note that Oc and TZc are irreducible. 

Theorem 13 (Waldhausen). — Let f : M ^ M' be a homotopy equivalence between ori- 
entable irreducible 3-manifolds with boundaries such that f takes the boundary of M onto the 
boundary of M' homeomorphically. Then f can be deformed to a homeomorphism M — )• M' 
by a homotopy which is fixed all the time on the boundary of M. (See |Mat03| . page 220, for 
a proof.) 

Applying Waldhausen's theorem, one obtains a homeomorphism ac from Oc to TZc which 
is identity on the sections in Sc- Gluing ac, one obtains a global homeomorphism form O to 
TZ. Moreover, according to Chazal and Cohen-Steiner's work |CCS05] (Corollary 3.1), since 
TZ and O are homeomorphic and TZ contains the medial axis of O, 7^ is isotopic to C □ 

Conclusion. — In this paper, we presented the first topological studies in shape reconstruc- 
tion from cross-sectional data. We showed that the generalization of the classical overlapping 
criterion to solve the correspondence problem between unorganized cross-sections, proposed 
by Liu et al. in J LBD"'"08] . preserves the homotopy type of the shape under some appropriate 
sampling conditions. In addition, we proved that in this case, the homotopy equivalence be- 
tween the reconstructed object and the original shape can be deformed to a homeomorphism. 
Even, more strongly, the two objects are isotopic. 
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A 

Homotopy Theory Basics 

In this section, we briefly recall some basic definitions that are used in the paper. 

Definition 12 (Homotopy). — A homotopy between two continuous functions / and g 
from a topological space A to a topological space Y is defined to be a continuous function 
i7 : A X [0, 1] — ^ y such that for all points x G A, H{x,0) = f{x) and H{x, 1) = g{x). f is 
said to be homotopic to g if there exists a homotopy between / and g. 

Definition 13 (Homotopy Equivalence). — Two topological spaces A and Y are homo- 
topy equivalent or of the sam,e hom,otopy type if there exist continuous maps / : A — )■ y and 
(7 : y — 7- A such that g o f is homotopic to the identity map idx and f o g is homotopic to 
idy- 

Definition I4 (Homotopy Groups, Fundamental Group). — Let A be a space with a 
base point xq G X. Let 5* denote the i-sphere for a given i > 1, in which we fixed a base 
point b. The i-dimensional homotopy group of A at the base point xq, denoted by 7ri(A, xq), 
is defined to be the set of homotopy classes of maps f : S^ ^ X that map the base point b to 
the base point xq. 

Thus if A is path-connected, the group 7rj(A, xq) is, up to isomorphism, independent of the 
choice of base point xq. In this case the notation 7rj(A, xq) is often abbreviated to 7rj(A). 

Let A be a path-connected space. The first homotopy group of A, 7ri(A), is called the 
fundam,ental group of A. 

Definition 15 (Simply-Connected Spaces). — The path-connected space A is called 
sim,ply- connected if its fundamental group is trivial. 

Definition 16 (Weak Homotopy Equivalence). — A map / : A — )■ y is called a weak 
homotopy equivalence if the group homomorphisms induced by / on the corresponding ho- 
motopy groups, /* : ■Ki{X) — )■ '7rj(y), for i > 0, are all isomorphisms. It is easy to see that 
any homotopy equivalence is a weak homotopy equivalence, but the inverse is not necessar- 
ily true. However, Whitehead's Theorem states that the inverse is true for maps between 
CW-complexes. 

Theorem I4 (Whitehead's Theorem). — If a map / : A — )■ y between connected CW- 
complexes induces isomorphisms /* : TTi{X) — )■ iTi{Y) for all i > 0, then / is a homotopy 
equivalence. 

Definition 11 ((Strong) Deformation Retract). — Let A be a subspace of Y . A ho- 
motopy i? : y X [0, 1] — 7- y is said to be a (strong) deformation retract of Y onto A if: 

— For all y G y, R{y, 0) = y and H{y, 1) G A. 

— For all X G A, ii{x, 1) = x. 

— (and for all x G A, il(x^ t) = x.) 

Definition 18 (Homeomorphism). — Two topological spaces A and y are /lomeomorp/iic 
if there exists a continuous and bijective map /i : A — )■ y such that /i^^ is continuous, the 
map h is called a homeomorphism from A to y. 
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Definition 19 (Isotopy). — Two topological spaces X and Y embedded in M are isotopic 
if there exists a continuous map i : [0, 1] x X — )■ M such that i(0, .) is the identity over X, 
i{l,X) = Y and for any t G [0, 1], i{t, .) is a homeomorphism from X to its image. The map 
i is called an isotopy from X to y. 

Let us recall the definition of the universal cover, and refer to classical books in topology 
for more details. 

Definition 20 (Universal Cover) . — Let X be a topological space. A covering space of X 
is a space C together with a continuous surjective map (j) : C ^ X such that for every x G X, 
there exists an open neighborhood U of x, such that (p~^{U) is a disjoint union of open sets 
in C, each of which is mapped homeomorphically onto U hy (f). A connected covering space is 
called a universal cover if it is simply connected. 

The universal cover exists and is unique up to homeomorphism. 

Lemma 4 (Lifting Property of the Universal Cover). — Let X be a (path-) connected 
topological space and X be its universal cover, and </> : X — t- X be the map given by the 
covering. Let Y be any simply connected space, and / : y — )■ X be a continuous map. 
Given two points x G X and y & Y with (f){x) = f{y), there exists a unique continuous map 
(7 : y — 7- X so that (p o g = f and (j){y) = x. This is called the lifting property of X. 

Since all the spheres Si of dimensions i > 2 are simply connected, we may easily deduce. 

Corollary 1. — For any (path-) connected space X with universal cover X, we have 7rj(X) = 
TTi{X) for all i > 2. 

We now recall the definition of the Hurewicz map hi : vrj(X) — )■ Hi{X). For an element 
[a] G TTi{X) presented by a : S"* — t- X, /ij([a]) is defined as the image of the fundamental class 
of 5* in Hi{S') under the map a* : Hi{S') -^ Hi{X), i.e., hi{[a]) = a,(l). 

Theorem 15 (Hurewicz Isomorphism Theorem). — The first non-trivial homotopy 
and homology groups of a simply-connected space occur in the same dimension and are 
isomorphic. In other words, for X simply connected, the Hurewicz map hi : ■Ki{X) — t- Hi{X) 
is an isomorphism for the first i with ttj (or equivalent ly Hi) non-trivial. 



